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XMM-NEWTON /SDSS: STAR FORMATION EFFICIENCY IN GALAXY CLUSTERS AND CONSTRAINTS ON THE 
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ABSTRACT 

It is believed that the global baryon content of clusters of galaxies is representative of the matter distribution 
of the universe, and can, therefore, be used to reliably determine the matter density parameter Q m . This 
assumption is challenged by the growing evidence from optical and X-ray observations that the total baryon 
mass fraction increases towards rich clusters. In this context, we investigate the dependence of stellar, and total 
baryon mass fractions as a function of mass. To do so, we used a subsample of nineteen clusters extracted 
from the X-ray flux limited sample HIFLUGCS that have available DR-7 Sloan Digital Sky Survey (SDSS) 
data. From the optical analysis we derived the stellar masses. Using XMM-Newton we derived the gas masses. 
Then, adopting a scaling relation we estimate the total masses. Adding the gas and the stellar mass fractions 
we obtain the total baryonic content that we find to increase with cluster mass, reaching 7-year Wilkinson 
Microwave Anisotropy Probe (WMAP-7) prediction for clusters with M500 = 1.6 x 1O 15 M . We observe a 
decrease of the stellar mass fraction (from 4.5% to ~1.0%) with increasing total mass where our findings for 
the stellar mass fraction agree with previous studies. This result suggests a difference in the number of stars 
formed per unit of halo mass, though with a large scatter for low-mass systems. That is, the efficiency of star 
formation varies on cluster scale that lower mass systems are likely to have higher star formation efficiencies. 
It follows immediately that the dependence of the stellar mass fraction on total mass results in an increase of 
the mass-to-light ratio from lower to higher mass systems. We also discuss the consequences of these results 
in the context of determining the cosmic matter density parameter Q m . 
Subject headings: cosmology: observations 



1. INTRODUCTION 

One of the key goals of observational cosmology is to de- 
termine the matter density of the universe (Q m ). One of the 
classic methods of inferring Q m is to use the lOortl d 1 9581) tech- 
nique adopting the mass-to-light ratios of galaxy clusters. In- 
dependently, one can assume that the ratio of baryonic-to-total 
mass of very massive galaxy clusters should closely match 
the ratio of the cosmo logical parameters, and thus Qh/ Qm ~ 
Mu/M, n , dWhite et al.lfl99l lEvrardl[T997l lEttori et al.ll200l 
lAllen et al.ll2008l) . 

In these two cases, it is widely believed that the global 
baryon content and the mass-to-light ratio of galaxy clusters 
are fairly representative of the matter distribution of the uni- 
verse, and can therefore be used, to reliably determine the 
cosmic matter density parameter. Contrary to expectations, 
this fundamental assumption is challenged by the growing 
evidence from optical an d X-ray observations that the to- 
tal baryon mass fraction dDavidet all 1 1 990t iLin et al.1 120031 : 
iGonzalez et alj I2007a|: | Giodini et alj I2009|) and the mass- 
to-light ratio (e.g. lAdami et all 11 998t~ Grrardi et al.l 120001: 



Bahca ll & Comerfordl 120021 iRines et al] 12004 iMuzzin et alj 
2001 ' increase towards rich clusters for a broad mass range of 



systems. Also the comparison between recent direct measure- 
ments of the baryon mass conte nt of galaxy cluster s and the 
prediction from the WMAP-7 dJarosik et alj 1201 11) presents 
controversial mass dependence of missing baryons (lAndreonl 
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l2Q"l0t iGiodini et al.ll2009T: IGonzalez et aO2007H) . 

Observational difficulties prevented recent studies from 
building a large sample of galaxy clusters where both masses 
and luminosities are homogeneously computed, what conse- 
quent ly affect the reliability of Q m estimates (iGirardi et alj 
120001) . From the optical point of view, the determination 
of cluster galaxy luminosities is limited by uncertainties re- 
lated to corrections for Galactic extinction and background 
galaxy contamination, extrapolation of the luminosity func- 
tion towards the faint-end, and the completeness of the sam- 
ple. Also, various methods can be applied to estimate clus- 
ter mass. The total masses can be inferred from either X-ray 
or optical data, under the assumption of hydrostatic equilib- 
rium. Estimates based on gravitational lensing do not require 
assumptions about the dynamical state of the cluster, but it 
encompasses substru ctures along the line of sight. Recent nu- 
merical simulations dMeneghetti et al.lfeOlOl) claim that com- 
bining weak and strong lensing data, the projected masses 
within r20(£] can be constrained with a precision of ~ 10%. 
However, deprojection of lensing masses increases the scat- 
ter around the true masses by more than a factor of two be- 
cause of cluster triaxiality. X-ray mass measurements have 
much smaller scatter (about a factor of two less than the lens- 
ing masses), but they are generally biased toward low values 
between 5% and 10%. 

The baryons in galaxy clusters consist of stars in clus- 
ter galaxies, intra-cluster light (IC L, stars that are not 
bounded to the cluster g alaxies; e.g.. IGonzalez et aTll2007ak 
Krick & Bernstein) 2007 ), and the hot intra-cluster gas (ICG, 



Sarazin & Bahcall 



19771) . The baryonic mass of a galaxy 



cluster is dominated by the ICG, the mass of which exceeds 

1 the radius at which the mean mass density is 200 times the critical density 
at the cluster redshift. 
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the mass of the former two components by a factor of ~ 6. 
Thus, to reliably compute this cosmological parameter from 
the baryonic-to-total mass ratio one should address all bary- 
onic components and not only the gas mass. In the past years, 
measuring individual baryon budget in galaxy clusters be- 
came possibJe_wkhj^ou^^ optical imaging data 
(e.g.. iDavid et al.lfl990l iLin et al.l 12003b iRoussel et alJ GOOOt 
iGonzalez et alj|2007bt lLagana et al.11200^ 7 The combination 
of the X-ray measured gas mass fraction (/ gas = M gas /M tot , ra- 
tio of the gas-to-total mass) with the optically measured stel- 
lar mass fraction (/* = M*/M tot , ratio of the stellar-to-total 
mass) can thus better constrain the total baryon mass fraction 
(/b — /* + /gas)- Notably, these observational measurements 
are essential to improve our understanding on the distributions 
of each of these components on cluster scales. Besides, stud- 
ies of structure formation require constraints on the baryonic 
components in galaxy clusters to understand star formation 
and metal enrichment processes. 

In this context, H IFLUGCS dReiprich & Bohringeri 12001 
iHudson et al.ll2010l) provides an X-ray flux-limited sample of 
64 nearby (z < 0.1) clusters selected from the ROSAT All- 
Sky Survey (RASS). Nineteen clusters in the HIFLUGCS 
sample have availab le Seventh Data Relea se of Sloan Digi- 
tal Sky Survey data dAbazaiian et al.ll2009l DR-7 SDSS;). It 
should be mentioned that the flux-limited samples might be 
biased towards more luminous clusters, leading to a larger 
fraction of cool-core clusters for lower mass systems and 
also they are not necess arily unbiased with respect to mor- 
phology dReiprichll2006l) . But the advantage of working with 
flux-limited samples is that the bias can be determined (e.g., 
Ilkebe et al. 112003 1 Stanek et"aT]|2006l: Eckert et al.ll20TTt) . 

The investigation of the total baryon mass fraction and 
mass-to-light ratio dependence on total cluster mass is one 
of the key goals to understand better the dynamical history 
of galaxy clusters and also the use of these systems to de- 
termine the cosmological parameter Q. m . Thus, in this paper, 
we analyze the X-ray XMM-Newton and optical DR-7 SDSS 
imaging data for the nineteen above mentioned clusters. It is 
important to stress that both luminosities and total masses are 
determined in a homogeneous way. Using the X-ray gas mass 
as a proxy to determine the total mass seems to be the more 
suitable technique since it is known that ga s mass may likely 
be unbiased mass proxy dOkabe et al.ll2010h . 

The paper is organized as follows. In Sect. [2] we describe 
the sample selection criteria providing X-ray and optical de- 
tails separately. We discuss our results in Sect. [3j and sum- 
marize our findings in Sect. [4] The scaling relation used 
to determine the total mass from the gas mass in given in 
Appendix [A] extra information about the clusters (such as 
the color-magnitude diagrams, CMD, and luminosity function 
fits) is given in Appendix iBl and a detailed discussion about 
the main systematic uncertainties involved to derive the stellar 
masses are presented in Appendix [Cl Throughout this paper 
we assume Q m = 0.3, £2a = 0.7 and, Hq = 70 km s _1 Mpc~'. 
Confidence intervals correspond to the 68% confidence level. 
We adopted these values to be c onsiste nt with the X-ray anal- 
ysis presented in IZhang et"aTl d201 lal) . Even so, it should 
be noted that the the WMAP-7 parameters for the standard 
ACDM model are Q m = 0.267 + 0.0288, O a = 0.734 + 0.029, 
and// = (71 ±2.5) km s" 1 Mpc" 1 . 

2. SAMPLE AND DATA ANALYSIS 

To constrain the mass dependence of the star forma- 
tion efficiency (SFE) and baryon mass fraction, it is im- 
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Fig. 1 . — XMM-Ne wton measured X-ray bolometric luminosity within the 
cluster radius f'500 in Zhang et al. (201 la) vs. redshift for the HIFLUGCS 
sample. The 19 clusters with SDSS imaging data are highlighted by circles. 



portant to use representativ e samples. The HIFLUGCS 
dReiprich & Bohringe r 2002) provides an X-ray flux-limited 
sample of 64 nearby (z < 0.1) clusters selected from the 
RASS. Nineteen clusters in the HIFLUGCS ha ve been ob- 
served in the DR-7 SDSS dAbazanan et alj|2009l) and are un- 
biased with respect to the original HIFLUGCS sample, having 
the same nature of flux selection, as can be seen in Fig.Q] 

2.1. X-ray data 

Detailed XMM-Newton data reduction can be found in 
IZhang et all (120091 1201 lal) including data screening, pointing 
source subtraction, background modeling, and spectral imag- 
ing analysis. 

For sample studies, one has to derive all quantities consis- 
tently within a characteristic cluster radius. We define them 
in terms of /-500, the radius at which the mean mass density 
within the cluster is 500 times the critical density at the clus- 
ter redshift. 

We can measure r^oo from the X-ray measured mass dis- 
tribution derived under the assumpt ion of hydrostatic equi- 
librium (HE; e.g. IZhang et al.l 120091) . However, the present 
sample contains relaxed an d non relaxed clus t ers wh ere the 
HE may not be valid (e.g., IZhang et al.ll2008l |20Tob . Thus, 
we used a scaling relation between the gas and the total mass. 
To do so, we u sed 41 dynamically relaxed groups and cluster s 
collected fromlVikhlinin et al.l d2006al) , [Arnaud et al.l d2007l) . 
iBohringer et al.l d2007l) . and iSun et ail ((20091) . We then fit a 
power-law to the data points and we used this fit to compute 
total masses (see Appendix[A]for more details): 



log 
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where E(z) = [On(l + z) 3 + (1 - On - 00(1 + zf + O0 1/2 , 
describing the redshift evolution of the Hubble parameter. 

We note that the following scaling relation is based on 
four different samples of groups and clusters what should re- 
duce the potential bias produced by sample selection. Be- 
sides, these 41 systems used to construct our scaling rela- 
tion are all relaxed systems, and thus the assumption of hy- 
drostatic equilibrium should be valid for these obejects and 
thus minimize the bias in total mass estimates for the clus- 
ters in our sample. In the way the total mass is computed, 
log M to t = A + B x log M gas , when the slope parameter (B in 
this case) is positive and lower than 1, the gas mass fraction 
encreases with total mass (see AppendixlAlfor more details). 

2.2. SDSS data 

We selected all HI FLUGCS clusters with DR-7 SDSS 
dAbazaiian et al.l 120091) data available. The DR-7 includes 
more accurate results than previous data releases, updated and 
more accurate magnitudes for all sources in the PHOTO cat- 
alog, and better accounting for sky background subtraction. 
Specifically, we use the dered table of magnitudes. Three 
clusters (Al 19, NGC507 and, A2063) are located on the bor- 
der of SDSS field of view, so that we had to exclude them. 
We ended up with 19 clusters spanning the redshift range of 
0.02 <z<0.1 listed in TableQ] 

We note that there are 2MASS dSkrutskie et al.ll2006h ob- 
servations for the HIFLUGCS clusters in which the K-band 
luminosity data have also been used to derive stellar masses. 
However, SDSS data allow us to make a cleaner selection of 
the cluster galaxies and to avoid fore- and background con- 
taminations, which are important to derive reliable luminosity 
functions (LF) and thus the stellar mass for our purpose. 

In order to obtain the stellar mass of galaxies we first con- 
structed color-magnitude diagrams (CMD; (g - i) x i) to 
statistically select cluster members within f-500. In this pa- 
per, we fit iteratively the red-sequence (RS), first identified in 
the 0.7 < (g - i) < 1.6 color interval and / < 18 magnitude 
limit. This interval was the one that best suited the majority 
of clusters in our sample. For four clusters (IIIZw54, A1367, 
MKW4, and A400) we had to extend the magnitude limit up 
to i < 16 to better define the slope of the RS. In this way, we 
select just the brighter and redder galaxies, what enables us to 
first identify the red-sequence on the CMD, fitting a linear re- 
lation to these galaxies. Then, to constrain the position and the 
slope of the RS on the CMD, we use the whole color and mag- 
nitude interval to re-do the linear fit until it converges. For red 
galaxies, we adopted ±0.3 mag width from the red-sequence 
best fit. Then, the blue galaxies were considered as the galax- 
ies that lay below the lower limit of the red-sequence (see Fig- 
ure lB.ll in AppendixIBl. To take into account fore/background 
contamination we computed the number of galaxies per 0.5 
bin of magnitude in an annular region beyond 8 x rsoo from 
the cluster center. We then fit the galaxy number counts in the 
20 < i < 16 range of magnitude using a power-law. The final 
cluster magnitude number counts were obtained by subtract- 
ing the background fit from the galaxy counts. 

For the selected blue and red members, w e co mputed 
the luminosity function in the /-band (see Fig. IB. II in Ap- 
pendix |B},gerformlrig_an analytical fit using two Schechter 
functions (ISchechterl 19761) with the usual parameters a (faint- 
end slope), O* (normalization) and M* (characteristic magni- 
tude). For all clusters in our sample but ZwC11215, A1650, 
and A2255 we used two Schechter functions to take into ac- 
count both bright- and faint-end contribution. With the im- 



provement of optical instruments, in the past years much work 
has been done to constrain the fai nt-population (-18 < M\ < 
-12) of nearby galaxy clusters (e.gJPennv & ConseUcell2008l : 
iRines & Gelleri 120081: iBoue et al] |2008j) Since the LF of 
our clusters showed a steep population faintward we fit 2 
Schechter functions to describe the overall behavior of galax- 
ies. The covariance between the characteristic magnitude and 
the power-law index is taken into account (see Appendix iBl 
for more details about LF fits). 

We considered jT-corr ection (K(z)) and evolutionary cor- 
rection (e(z)) supplied by iPoggiantil (1 1 9971) to correct M* val- 
ues. To derive the total luminosity we integrated the lumi- 
nosity function assuming a lower luminosity limit. We thus, 
assumed the corresponding luminosity of Mj = -14.0 for all 
clusters, in order to integrate the luminosity function down 
to the same flux limit. To obtain the stellar masses, we 
adopted differe nt mass-to-light ratio fo r ellipticals and spi- 
rals, ta ken from lKauffmann et al.l (120031) assuming a lSalpeteri 
(119551) initial mass funct ion (IMF), as stated in Equation 21 
from lLagana et alj (120081) . The stellar mass estimate is tied to 
the choice of IMF. Changing the IM F scales the ste llar mass 
estimate by a fi xed factor, e.g., from a lKroupal (120011) IMF to a 
ISalpeted (119551) IMF with a cut-off at 0. 1M results in a facto r 
of 2 increase in the stellar mass (e.g..,Kauff mann et al.ll2003l) . 
These procedures were also applied in lZhang et al.l d201 lbl) . 

When there is a systematic offset between the spectroscopic 
and photometric magnitu des, the stellar mass e stimate can 
be systema tically biased dKauffmann et alJ l2003). The cal- 
ibration in iFritz et all (1201 ll) shows that the stellar masses 
computed from fiber-aperture magnitudes using DR-7 SDSS 
photometric data are lower by ~ 0.15 dex than the values 
computed from spectroscopic data. Our sample is in simi- 
lar redshift and mass ranges as their sample. We, therefore, 
corrected this bias in our stellar mass estimates using best- 
fit relation between stellar masses computed from the DR-7 
SDSS photometric data an d from the spectro scopic data (see 
the right panel in Fig. 4 in IFritz et alj|201 ll) . All values for 
stellar masses (M™ 1T ) and optical luminosities C^SSD given in 
Table [2] are already corrected and are derived on the basis of 
7"500- We also give the a and M* parameters for the Schechter 
fit in this table. A detailed discussion about the impact of 
the main systematic effects on M c ° n and L z °^ are presented in 
AppendixICl 

3. RESULTS: FROM STELLAR MASS FRACTION TO 
THE COSMOLOGICAL MATTER DENSITY 
PARAMETER 

To fit the trend in our results, we apply the BCES regres- 
sion fitting metho d taking into account mea surement errors 
in both variables (Akritas & Bershady 1996). Through this 
paper, we apply the BCES regression fitting taking into ac- 
count me^suremejit_OTO£s in_both variables and their covari- 
ance (lAkritas & Bershadylll996l) . Confidence intervals corre- 
spond to 68% confidence level. 

Figure [2] shows the stellar, and total baryon fraction as a 
function of total m ass. We also compared our result to some 
previous analyses (iLin et al.ll2003t [Giodini et alJl2009l) . We 
discuss these results in detail in Section lXTl 

In Figure[3] we show the total mass-to-optical light ratio as 
a function of total mass. This result is specifically discussed 
in Section USl In Section l3~3l we discuss the consequences of 
determining the matter-density parameter using galaxy clus- 
ters. 
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TABLE 1 

HIFLUGCS clusters with DR-7 SDSS available data 



Name 


K.A. 




reasnilt 


rsao (Mpc) 


A/f ha 

M g;ls (10 


13 A A \ 


a a / 1 r 




A85 


10.46 


-9.30 


0.0556 


1.216 


± 0.058 


7.46 ± 


0.34 


5.68 ± 


0.81 


A400 


44.42 


6.02 


0.0240 


0.712 


± 0.034 


1.00 ± 


0.04 


1.07 ± 


0.15 


lllZw54 


55.32 


15.40 


0.031 1 


0.731 


± 0.035 


1.10 + 


0.21 


1.18 ± 


0.17 


A 1367 


176.01 


19.95 


0.0216 


0.893 


± 0.043 


2.24 + 


0.07 


2.11 ± 


0.30 


MKW4 


181.11 


1.90 


0.0200 


0.580 


± 0.028 


0.47 + 


0.02 


0.58 ± 


0.08 


ZwC11215 


184.42 


3.66 


0.0750 


1.098 


± 0.052 


5.40 + 


0.29 


4.34 ± 


0.62 


A 1650 


194.67 


-1.76 


0.0845 


1.087 


± 0.052 


5.36 ± 


0.90 


4.28 ± 


0.61 


Coma 


194.90 


27.96 


0.0232 


1.278 


± 0.061 


8.18 + 


0.62 


6.21 ± 


0.89 


A1795 


207.22 


26.59 


0.0616 


1.118 


± 0.053 


5.57 ± 


0.14 


4.46 ± 


0.63 


MKW8 


220.18 


3.46 


0.0270 


0.715 


± 0.034 


1.02 ± 


0.14 


1.10 ± 


0.16 


A2029 


221.13 


5.74 


0.0767 


1.275 


± 0.061 


9.17 + 


0.33 


6.81 ± 


0.97 


A2052 


229.19 


7.02 


0.0348 


0.875 


± 0.042 


2.13 + 


0.11 


2.03 ± 


0.29 


MKW3S 


230.47 


7.71 


0.0450 


0.905 


± 0.043 


2.44 ± 


0.10 


2.29 ± 


0.33 


A2065 


230.60 


27.71 


0.0721 


1.008 


± 0.048 


4.12 ± 


1.29 


3.35 ± 


0.48 


A2142 


239.58 


27.23 


0.0899 


1.449 


± 0.069 


15.10 ± 


0.88 


10.26 ± 


1.47 


A2147 


240.57 


15.97 


0.0351 


1.064 


± 0.050 


4.35 + 


0.43 


3.63 ± 


0.52 


A2199 


247.16 


39.55 


0.0302 


0.957 


± 0.045 


2.90 ± 


0.23 


2.64 ± 


0.38 


A2255 


258.12 


64.07 


0.0800 


1.072 


+ 0.051 


5.13 + 


0.20 


4.08 ± 


0.58 


A2589 


350.99 


16.78 


0.0416 


0.848 


± 0.040 


1.93 ± 


0.12 


1.88 ± 


0.27 



TABLE 2 
Photometric results 



Name 




Brit 

«i 


;ht-end 




Faint-end 


m; 




M cor,- (10 i2 Mo) 


/* 




A 


A85 


-1 


42 ± 0.09 


-22.65 ± 0.32 


-1 


.78 ± 0.080 


-20 


.05 ± 0.51 


2.32 ± 0.24 


6.93 ±0.71 


0.0122 ± 0.0022 


0.144 ±0.020 


A400 


-1 


.15 ±0.06 


-22.16 ±0.19 


-1 


.77 ± 0.002 


- 18 


.72 ±0.31 


1.61 ±0.20 


4.93 ± 0.62 


0.0459 ± 0.0087 


0.139 ±0.016 


IIIZw54 


-1 


.19 ±0.08 


-21.82 ±0.65 


-1 


.71 ±0.044 


- 18 


.66 ± 0.64 


1.38 ±0.18 


3.96 ± 0.52 


0.0337 ± 0.0065 


0.127 ±0.023 


A1367 


-1 


.27 ± 0.06 


-22.66 ± 0.58 


-1 


.84 ± 0.020 


- 18 


.18 ±0.31 


1.59 ±0.27 


4.75 ± 0.80 


0.0225 ± 0.0050 


0.129 ±0.016 


MKW4 


-1 


.10 ±0.16 


-22.00 ± 0.36 


-1 


.78 ± 0.020 


- 18 


.01 ±0.13 


0.31 ±0.03 


0.93 ±0.11 


0.0162 ± 0.0030 


0.097 ± 0.012 


ZwC11215 


-1 


.64 ± 0.03 


-23.13 ±0.31 










2.71 ±0.33 


8.14 ± 1.00 


0.0187 ±0.0035 


0. 


143 ± 0.019 


A 1650 


-1 


.67 ± 0.01 


-23.25 ±0.31 










2.24 ± 0.38 


6.67 ± 1.12 


0.0156 ± 0.0034 


0. 


141 ± 0.028 


Coma 


-1 


.26 ± 0.03 


-21.95 ±0.30 


-1 


.98 ±0.011 


- 18 


.06 ± 0.08 


4.76 ± 0.54 


14.00 ± 1.59 


0.0226 ± 0.0041 


0. 


154 ± 0.022 


A1795 


-1 


.22 ±0.01 


-21.41 ± 0.13 


-1 


.46 ± 0.008 


- 19 


.34 ± 0.41 


2.27 ± 0.34 


6.99 ± 1.04 


0.0157 ± 0.0032 


0. 


141 ±0.018 


MKW8 


-1 


.30 ±0.11 


-21.79 ±0.58 


-1 


.74 ± 0.009 


- 18 


.68 ± 0.70 


0.61 ± 0.09 


1.76 ±0.26 


0.0161 ± 0.0033 


0. 


109 ± 0.019 


A2029 


-1 


.17 ± 0.07 


-21.69 ±0.17 


-1 


.66 ± 0.020 


- 19 


.84 ±0.12 


3.18 ±0.41 


9.49 ± 1.22 


0.0139 ±0.0027 


0. 


148 ± 0.020 


A2052 


-1 


.27 ± 0.10 


-22.10 ±0.48 


-1 


.70 ±0.011 


- 18 


.59 ± 0.09 


1.37 ±0.14 


4.16 ±0.42 


0.0205 ± 0.0036 


0. 


125 ±0.016 


MKW3S 


-1 


.36 ± 0.03 


-22.14 ±0.48 


-1 


.57 ± 0.041 


- 19 


.12 ±0.25 


3.13 ±0.19 


9.39 ± 0.80 


0.0410 ± 0.0068 


0. 


147 ± 0.017 


A2065 


-1 


.10 ±0.02 


-22.00 ±0.11 


-1 


.62 ± 0.027 


- 20 


.08 ±0.12 


2.55 ± 0.27 


7.49 ± 0.57 


0.0224 ± 0.0036 


0. 


145 ± 0.042 


A2142 


-1 


.00 ± 0.02 


-21.84 ±0.86 


-1 


.68 ± 0.024 


- 20 


.56 ± 0.11 


3.51 ±0.19 


11.00 ±0.96 


0.0107 ± 0.0018 


0. 


158 ± 0.023 


A2147 


-1 


.42 ± 0.05 


-21.93 ±0.12 


-1 


.76 ± 0.015 


- 18 


.93 ±0.12 


2.24 ±0.31 


6.98 ± 0.85 


0.0192 ± 0.0036 


0. 


139 ±0.021 


A2199 


-1 


.18 ±0.03 


-21.80 ± 0.20 


-1 


.86 ± 0.009 


- 18 


.28 ±0.56 


1.51 ±0.27 


4.72 ± 0.56 


0.0179 ± 0.0033 


0. 


128 ±0.018 


A2255 


-1 


.45 ± 0.03 


-23.10 ±0.18 










2.00 ± 0.26 


6.18 ±0.79 


0.0151 ±0.0029 


0. 


141 ±0.019 


A2589 


-1 


.09 ± 0.02 


-21.65 ± 0.20 


-1 


.58 ± 0.010 


- 18 


40 ± 0.86 


1.36 ±0.14 


4.19 ±0.43 


0.0223 ± 0.0039 


0. 


125 ±0.017 



3.1. Star formation efficiency 

As clusters of galaxies are the largest virialized systems in 
the Universe, the baryon budget detected in clusters of galax- 
ies should be relatively representative of the baryon mass frac- 
tion estimate of the Universe, when cons idering very massive 
halos (e.g.. lEttori et al.ll2003b IXilen et al.ll2008h . 

Thus, the total baryon fraction, ft,, would simply be given 
by Mb = /b M tot , where Mb is the sum of M gas + M*, and 
M tot is the hydrostatic mass. Using our sample of 19 clusters, 
we investigated the dependence of /*, and f\, as a function of 
cluster mass. 

In Figure |2] we observe a decrease of the stellar mass frac- 
tion, and an opposite trend, although less sharply for the gas 



mass fraction. More massive clusters have gas mass fraction 
two times higher and stellar mass fraction around four times 
lower than less massive clusters. We fit the behavior of the 
stellar and total baryon-mass fractions with mass as power 
laws. Within r^m the power-law fit for the stellar-mass frac- 
tion is /* = 10 ( - L54±01 °) x [M 5O o/lO 14 M ] ( - a36±o - 17) , and the 
total baryon-mass fraction is given by f\, = io ( ~°- 930±0018 ) x 
[M 5O o/lO 14 M ] (O136±OO28) . At this point, we call the attention 
to the fact that the observed increase of the ft, with total mass 
can be partially due to the assumption that the gas mass in- 
creases with total mass as discussed in AppendixlAl 

Our study shows a decrease of the stellar mass fraction 
from 4.5% to approximately 1.0%, and taking into account 
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the errors involved, one should notice this dependence devi- 
ates from a fiat distribution within ~95% confidence. This re- 
sult suggests that, at least on cluster scales, the number of stars 
formed per unit of halo mass between low and high-mass clus- 
ters are different, that is the star formation efficiency (SFE) de- 
pen ds on the environment (as already addressed in several pa- 
pers IGonzalez et al.l2000tlLin et al.l200 3t lLagana et al.l2008b 
lEttori et al.l 120091: IZhang et alj|201 lbl) . In massive clusters, 
more hot gas and dark matter is settled in the deep potential 
wells than the amount of individual galaxies. In low-mass 
systems, the accretion of individual less massive galaxies is 
more important and more low entropy gas is brought in to 
form stars. Both processes result in lower star formation effi- 
ciency for more massive systems. 

As a consequence of the dependence of the baryon mass 
fractions on total mass, one does not observe a "fiat" /b dis- 
tribution as a function of the total mass on cluster scales. As 
shown in Figure|2] the total baryon mass fraction is lower than 
the WMAP-7 predic t ed cosmic fraction, in line with previous 
studies, e.g.. lEttoril (12003b : iMcCarfhv etail (12007b . but see 
ISimionescu et al.l (1201 lb . In particular, the discrepancy be- 
tween the observed and the fit value becomes larger with de- 
creasing mass reaching ~ 5.5cr for systems with masses below 
3 x 10 14 M o , and ~ 3<x for systems with their masses higher 
3 x 1O 14 M . The disagreement between o bservations and 
WMAP result has been already pointe d out: IGonzalez et all 



(2007)3) presented a 3.2<x discrepancy. IGiodini et al 
found a 5cr disparity, and more recently lAndreon 



(2009) 
(2010) 



showed a 6cr distinction. These studies analyzed cluster and 
groups sample in a similar mass range as the one used in this 
work (5xl0 13 M^-10 15 M,A IGiodini et ap(|2009j) an alyzed 91 
object s spaning a range in M500 of ~ 10 13 - 10 15 M^. [Andreonl 
(2010) investigated 52 clusters and groups in the same mass 
range as the latter authors, and IGonzalez et alj d2007bb con- 
sidered 24 galaxy clusters in a broader interval of mass, rang- 
ing from M500 = 10 12 M Q up to M 5() o ~ 8 x 1O 14 M . 

It is worth noting that the most massive cluster in our sam- 
ple has a total mass of M ~ 1O 15 M and its baryon mass 
fraction already reaches ~96% of the cosmic value. Extrapo- 
lating our fit for the total baryon mass fraction, fb reaches the 
WMAP-7 prediction for clusters with M > 1.6 x 10 15 M o . We 
show the ratio between observed and WMAP-7 values for fc> 
as a function of total mass in the lower panel of Fig. [2] Hydro- 
dynamical simulations claim that baryons in clusters have un- 
dergone a large depletion during the formation of these struc- 
tures. 

Although lWhite et aD (1 1993b claim that within a sufficiently 
large radius the mean ba ryon fraction must take the global 
value, iFrenk et ail d 1999b assert that the depletion parameter, 
that is the ratio between the observed and the predicted baryon 
mass fraction is T - 0.925, which is in line with mo re re- 
cebt estimate (e.g.. lEttori et al.l2006tlYoung et al.l201 lb . This 
means that even considering a large radius, e.g. ^200, some 
baryons are missing also in massive clusters. 

The comparison of the observed baryon b udget to the 
WMAP-7 cosmic fraction (Jaros ik et al.l 1201 lb consistently 
suggests that the amount of mis sing baryons inc reases to- 
ward low-mass systems. Recently. iDai et all (12010b proposed 
that the baryon loss mechanism is primarily controlled by the 
depth of the potential well of the system. For deep poten- 
tial wells, such as rich clusters, baryon loss is not signifi- 
cant, however, for lower-mass clusters, baryon loss becomes 
increasingly important. Such a mechanism could be the pre- 



heating of baryons before they collapse. In this picture, the 
baryons never fall into groups and low-mass clusters, but re- 
main well beyond r2oo- Alternatively, the gas falls into the 
potential-wells and is subsequently removed by feedback pro- 
vided from SNe and AGNs, which should also play a major 
role for low-mass systems. 

Some missing baryons are suggested, which can partially 
account for the discrepancy between the observed and the pre- 
dicted cosmic baryon mass fraction. Intra-cluster light (ICL) 
is one of the most important missing baryons. Observational 
results have shown that the ICL can acco unt for 6%-22% of 
the to tal cluster light in the r-band (e.g.. iRrick & Bernstelnl 
120071: IGonzalez et al.l2007bl:lPierini et al.l2008b . and being an 
order of magnitude lower, it would account for ~ 2% in the 
total baryo n mass fraction. Recen t numerical simulations per- 
formed by iPuchwein et al.l (12010b showed that the amount of 
gas and ICL removed by active galactic nucleus (AGN) heat- 
ing from the central regions of clusters and driven outwards 
(r > rsoo) depends on clus t er mas s, being higher in low-mass 
systems. IPuchwein et al.l (12010b estimate that for a cluster 
with total mass M - 4 x 10 13 M o the ICL account for 58% 
of the total stars. Considering this latter result, there is no sig- 
nifican t missing baryons as suggested by e.g. lAfshordi et al.l 
(12007b . Just the ICL and some part of the ICG would be 
sufficient to explain the difference between observations and 
cosmological prediction. However, such high amount of ICL 
would be visible in current observations. 

Another reasonable possibility is a dependence 
of the baryon bud g et as a function of radius (e.g., 
iHenriksen & Mamonl 119941: lEttori & Fabian! 11999b . It 
might be that in the limiting radius of X-ray observations 
(rsoo) t he total baryon bu dget is still underestimated. Al- 
though |Allenet al. (2002) show that the gas mass fractions 
in the clusters asymptote towards an approximately constant 
value at a radius /-2500, most of the numerical simulation 
results and observational data show a n evident increa s e of the 
baryonic content with cluster radius. Vikhlinin et Id] (l2006bb 
showed that the gas mass fraction i ncreases with radius as a 
power law of overdensity. Recently, ISimionescu et al.l (1201 lb 
showed that the baryon mass fraction increases dramatically 
with cluster radius for the Perseus cluster. This cluster is 
one of the few clusters having deepest X-ray data to measure 
the baryon mass fraction up to the virial radii. They found 
that, for the Perseus cluster, the gas mass fraction alone 
(without considering the stars) reaches the baryonic cosmic 
mean value at about half of r2oo and at r2oo it is ~ 38% above 
the mean baryonic value. However, it is important to have 
in mind that this is an isolate result and it is important to 
analyze the dependence of / gas (determined within the virial 
radius) on total mass for a sample of clusters. For instance, 
iReiprichl (|2003) studied the gas mass fraction within ^200 for 
106 clusters, and Perseus is among the five clusters with the 
highest / gas . This study can be readdressed using Suzaku 
data. 

3.2. Total mass-to-optical light ratio 

In this work we used independent X-ray and optical meth- 
ods to compute the total mass-to-optical light ratio for our 
sample, and both optical luminosities and total masses were 
homogeneously determined. Compared to other estimates, 
the determination of total masses using X-ray is a robust 
and a preferable method since it involves lower uncertainties 
dNagai et al.ll2007llMeneghetti et aLlHoTol) . 

For ACDM cosmological model, the mass-to-light ratio 
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Fig. 2. — Upper panel.ToXal baryon-mass fraction (black circles), gas-mass fraction (green triangles) and stellar-mass fraction (blue stars) as functions of the 
total mass. The black and the blue solid lines represent power-law fits for the total baryon and stellar mass fractions as a function of total mass, respectively. The 
horizontal lines on the top of the plot represent the WMAP-7 result (/WMAP-7 = 0.171 ± 0.009 IJarosik et aljj201 1|) with 1 sigma error. For comparison, we also 
show the best-fit for the baryon-mass fraction from Lin et al. (2003) and for the stellar mass fraction from Giodini et al. (2009). Lower panel: Ration between 
observational and WMAP-7 baryon mass fraction as a function of total mass. 



(M/L) in high-mass systems is supposed to be approximately 
independent of halo mass. It has been generally found that 
M/L increases with halo mass but there is an approximate 
plateau in M/L v alues for the richest bound systems (e.g ., 
iDavid et al.lfl99l iCirimele et al1[l997l LSheldon et al.ll200l . 
The existence of this plateau is sometimes taken as evidence 
that the measured values of M/L do indeed represent the uni- 
versal value. 

However, recently, most authors have found a dependence 
between the M/L and the total mass, even for massive clus- 
ters of galaxies. Assuming a power-law relation, M500/L+ °c 
M 5 q , those authors have found a to be in the range of 0.2- 
0.4, in both optical and near- i nfrared bands, over a large 
mass range (lAdami et al.lll998b iBahcall & Comerfordl 120021: 
iRines et al.ll2004l) . 

In Figure [3] we show the dependence of the total mass-to- 
optical light ratio (Msoo/i*) on total mass. The best-fit for the 
mass-to-optical light ratio follows M 50 o/i* = 10 (2 02±0 I0) x 
[M 50 o/10 14 M ] (0 - 361±0169) . Considering the sample analysed 
in this work, a plateau is not clearly present for the massive 
end, and the best-fit leads to a Msoo/i* = 241M O /L for clus- 
ters with total mass of 1O 15 M . 

We see from Figure [3]that non cool-core clusters present a 
higher dispersion when compared to cool-core clusters. This 
may indicates that the dynamical state of the cluster should be 
considered when analyzing the mass-to-light ratio, because 
it might follow that for a large and homogeneous sample of 
cool-core clusters a plateau becomes more evident. Thus, in 
order to better investigate whether there is a plateau at the 
massive end, further studies should consider a complete flux- 
limited sample of cool-core clusters. 

The dependence of M500/L* shown in Figure [3] with clus- 
ter mass is a direct consequence of the varying star formation 
efficiency in galaxy clusters. Since low-mass systems have 



O NCC 




Fig. 3. — Total mass-to-light ratio as a function of total mass of the system. 
The black, red and green colors denote strong, weak cool-core and non cool- 
core clusters, respectively. 



higher star formation efficiency compared to a massive clus- 
ter, the optical light will not vary uniformly and thus, the total 
mass-to-optical light shows a dependence on cluster mass. 

3.3. The cosmological parameter Q m 

The results presented in the Sections [3.11 and [3~2l have im- 
portant implic ations on the cosmological use of galaxy clus- 
ters. Using the lOortl dl958l) technique, one can determine the 
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matter density parameter following Q m = (M/L) x Pl/Pc, 
where p c is the critical density, and p^ is the typical lumi- 
nosity density of the universe. A fundamental assumption 
in this application is that M/L of clusters is a fairly rep- 
resentative measurement of the universal value. However, 
the mass-to-light ratio has been found to increase with clus- 
ter mass. Assuming the mass-to-light ratio range shown in 
Figure [3] (from 70 M /L up to 300 M /L ), and adopting 
ji ~ 2.12 x 1 8 /iion Lp/Mpc for the luminosity density (see 
Table 3 from iBlanton et all 120031) . we have the total matter 
density parameter spanning the range 0.07 < Q m < 0.27. 
Taking into account the errors involved in this method to de- 
rive both total mass (~ 10%) and the optical light (a mean 
value of about 20%), we can determine f2 m within ~ 14% of 
uncertainty. It should be noted that since we did not reach 
the claimed plateau in mass-to-light ratio towards the mas- 
sive end, the upper limit obtained here for the matter density 
parameter (Q m = 0.3), although consist ent with the WMAP - 
7 concordance value within the errors dJarosik et al.l 1201 lb . 
could still be an underestimate of Q m . 

Another way to determine the matter density parameter is 
considering that the ratio of baryonic-to-total mass in clusters 
should closely match the cosmological param eters Qb/Q m . 
Thus adopting Cl b h 2 = 0.02260 + 0.00053 dJarosik et al.l 
|201 ll) . and the total baryon fraction computed in this work 
(and listed in Tabled we obtained 0.15 < O m < 0.30, a nar- 
rower interval compared to the previous method. This shows 
that the stellar and gas mass estimates from optical and X- 
ray data provide better determination of Q m when compared 
to the use of M500/L*. The first advantage is that to derive 
the gas mass one do not assume hydrostatic equilibrium, and 
thus the dynamical state of the cluster will not bias this deter- 
mination. Then, the unce rtainties involved i n this determina- 
tion a re around 5%-10% dNagai et al.ll2007tlMeneghetti etail 
120101) . Also, X-ray gives a 3D view of the cluster when de- 
projection is done, while total light is always derived from a 
2D projection of galaxies in the CMD. There are also the dif- 
fic ulties already pointe d out in the Section Q] and discussed 
in iGirardi et al.l d2000l) related to the luminosity estimates. 
Moreover, there has been some concern with the Mj L method 
because recent numerical simulations suggest that light is a 
biased tracer of dark matter, providing a biased val ue for D. m 
dOstriker et al.l2003l) . The Oort method|OorJ dl958l) necessar- 
ily focus on denser regions of the galaxy distribution, and thus 
the estimate of O m relies on the assumption that the galaxy 
population in these regions is representative of the universe 
as a whole. However, it is well-known that the galaxy pop- 
ulation in clusters differs from the field population since the 
work from|DressleJ (119801) . Since rich clusters of galaxies en- 
compasses large regions (of several Mpc), observations that 
extend up to the virial radius (r ~ r2oo) should, in principle, 
contain a sufficient collapsed region and provide a fair mass- 
to-light ratio. Nevertheless, it is important to stress that in this 
work we just examine a cluster volume of radius r^oo, what 
can also affect our estimates. 

It is worth mentioning the effect of the ICL on both deter- 
minations. In both cases, th e ICL would contribute to lower 
the Q m estimates. Using the lOorj (119581) technique, the ICL 
would account for the total light and then this ratio would 
be smaller than the ones obtained here. Thus, if on the one 
hand the upper limit for Q. m could still be higher when reach- 
ing the plateau for massive clusters; on the other hand, the 
inclusion of the ICL can contribute to a lower estimate of 
the matter density parameter. Furthermore, using the total 



baryon-to-mass fraction, the matter density parameter is given 
by O m ~ fib x M tot /Mb, and the ICL would account for the 
total baryon estimates (Mb). Consequently one would obtain 
a lower determination of the matter density parameter. Thus, 
in both cases, the derived values in this work are biased up for 
not considering the ICL contribution. 

4. CONCLUSIONS 

We studied a representative sample of 19 clusters, and in- 
vestigated this sample as a robust local reference to constrain 
the baryon mass fraction and its mass dependence on cluster 
scales for future studies. Our conclusions are as follows. 

• We obtained optical luminosities and total masses for 
a sample of 19 clusters in a homogeneous way. Thus, 
we computed total mass-to-optical light ratio using two 
independent methods, what present a considerable ad- 
vantage over previous studies. 

• The total mass-to-optical light ratio decreases toward 
low-mass systems, following the relation M500/L* = 
lO (2 O2±o lO) x[M 5 oo/lO 14 M ]( a361±0169 ). This result is a 
direct consequence of the varying SFE in cluster scales. 
In this work the Msoo/i* varies from ~60 M /L up 
to almost 300 M /L . Total mass-to-optical light ratio 
does not show evidence of a flattening, and the best-fit 
leads to a M500/-L* = 241M /L for clusters with total 
mass of 1O 15 M . Our results indicate that this flattening 
towards massive clusters should be more evident when 
considering a flux-limited sample of cool-core clusters. 
We observed that non cool-core clusters present a large 
dispersion when compared to more relaxed systems and 
may introduce some bias in the final dependence. 

• Both /b and /* show a mass dependence. Within 
r 50 o, the power-law fits are /+ = io<- L54±OIO > x 
[M 5O o/lO 14 M ]<- a359±ai7O >, and /„ = i (- <»±o.oi8) x 
[M 5O o/lO 14 M ] (OI36±OO28) for the stellar, and total 
baryon mass fraction, respectively. The observed in- 
crease of the /b with total mass can be partially due 
to the increase of the gas mass with total mass as dis- 
cussed in Appendix lAl Toward the high-mass end, the 
value gradually approaches the WMAP-7 prediction, 
reaching it at M500 = 1.6 x 1O 15 M , when extrapolating 
our fit for the total baryon mass fraction. However, the 
most massive cluster in our sample has a total mass of 
M500 = 1 -026 x 1O 15 M and its baryon mass fraction al- 
ready reaches ~96% of the cosmic value. The amount 
of missing baryons with respect to the WMAP-7 pre- 
dicted fraction increases toward lower systems. This 
can be compensated by ICL and a significant part of 
the gas that is driven outwards (r > rsm) due to AGN 
feedback, especially in lower mass-systems. 

• The SFE is lower in more massive clusters, increas- 
ing towards the low-mass end. The rapid decrease ob- 
served in the stellar mass fraction as a function of total 
mass suggests a significant change in the efficiency of 
star formation wi th cluster mass (as already addressed 
in several papers: iGonzalez et al.ll2000t OneFaIll200l 
lLagana et alj|2008t jEttori et al.1120091 

• We derived the matter density parameter using the 
lOortld 19581) technique and also from the bary on-to-total 
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mass ratio. Using these two approaches, we obtained 
0.07 < Q m < 0.3 and 0.15 < Q m < 0.27, respec- 
tively. Using the baryon-to-total mass ratio to compute 
Q m seems to give narrower range and more accurate 
values. Since we did not consider the ICL contribution, 
the obtained values for the matter-density parameter are 
possibly biased low. 
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APPENDIX 
A. M GAS -M m RELATION 

Since total mass estimates are derived under the assumption of hydrostatic equilibrium (EQ), and the present sample contains 
clusters in a wide variety of dynamical states (here separated in strong cool-core, weak cool-core and non cool-core clusters), 
the EQ assumption may not be valid in all cases. Thus, we preferred to derive total masses using the scaling relation. Recent 
observ ational investigations using a variety of cluster samples have demonstrated that the gas mass is indeed a low-scatter mass 
proxy (lMaughan1l2007l : lAxnaud et"ai1l2007l). To construct a scaling - relation between tot a l and gas mass, as show n in Fig.Q] we 
used 4 1 dynamically relaxed clusters from I Vikhlinin e t al. ( 2006ah. [Arnaud et alJ (120071) . iBohringer et all (120071) . and lSun et alJ 
(120091) . where the dynamical equilibrium assumption can be fairly applied to derive total mass. Then, using M tot - M gas relation 
stated in Eq.Q] we can recove r the total mass of our sample. In this way, for the non-relaxed clusters, the total mass can be quite 
accur ately determined (~10%, Nagai et al. 2007; Zhang et al. 2008). Also, according to numerical simulations (e.g., |Nagai et alJ 
120071) . the total ICM mass is measured quite accurately (<6%) in all clusters, irrespective to the dynamical state. 

We also should highlight that the way the total mass was computed from the scaling relation, log M tot = A + B x logM gas , one 
obtaines M tot = 10 A X M^ as , and / gas = M gas /M tot = IQ~ A ' B x (M tot ) (l - B)/B . Thus, being / gas oc {M lot ) {l - B>/B , when B is positive 
and lower than 1, the gas mass fraction increases with total mass. Since we obtained B = 0.827, we assume / gas oc M[' ot 2 . 




Fig. 1. — Total mass as a function of gas mass for 41 groups and clusters. The continuous line is the power-law fit for the data. 
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B. COLOR-MAGNITUDE DIAGRAMS AND THE LUMINOSITY FUNCTION FITS 

The el liptical galaxies in individual clusters form a red-sequence (RS) with a well-defined slope and small scatter dBower et al.l 
Il992bla[) . Specifically, the existence of the red-sequence at higher redshifts indicates that cluster elliptical galaxies are a passively 
evolvi ng population in which the reddening of massive galax ies is the result of a mass-metallicity relation rather than an age 
effect (Kodama & Arimotoll 19971: iKauffmann & Charlotlll998h . From the obse rvational point-of-view, studies have sho wn that 
for nearby clusters (z < 0.05) the slope of the red-sequence is almost constant (Gladde rs et al.|[l998t iRomeo et al.ll2008l) . Since 
we have almost half of the clusters in our samp le wi th redshifts higher than the above mentioned redshift, we did not fix the slope 
of the RS but left it to vary. We show in Figure IBTTI the red-sequence for all clusters in our sample. Even not fixing the slope, we 
obtained well-defined red-sequence in which the mean slope is b — -0.036 + 0.010. 

For all galaxies belonging to the cluster, we build the luminosity function. The luminosity function is the distrib ution of all 
morp hological types of cluster galaxies over the magnitude. This distribution is generally fit by the Schechter function dSchechterl 
119761) . In Figure IB. II we present the luminosity function fits for all nineteen clusters. In particular, the slope of the faint end 
of the LF is a direct indicator of the importance of dwarf galaxies, which are expected to be more fragile in the environment of 
clusters. The great majority of studies of th e LF indicate faint-end slopes in the range of -0.9 to -1.5, but these mostly did not 
reach very faint magnitudes (see Table 1 in iDe Propris et al1l2003L and references therein). Recent studie s using deep ima ging 
has shown deep estimations of the faint-end slope within the range -2.29 < a < -1.07 (see Table Al from lBoue et al.1120081 and 
references therein). Our mean values for the slope of the bright population and of the faint-end are < a\ >= -1.28 + 0.18 and 
< »2 >= -1.72 + 0.13, respectively. For the characteristic magnitude we obtained M\ - -22.16 + 0.50 and M* 2 = -18.97 + 0.85 
for the bright- and for the faint-end, respectively. 

C. IMPACT OF SYSTEMATIC UNCERTAINTIES IN THE PHOTOMETRIC RESULTS 

One of the basic observational re sult of the present study is the stellar masses estimate. To correct determine the stellar mass 
we followed the description in Sec. 12.21 and here we will discuss the principal steps that account for the systematic uncertainties 
involved in the total optical luminosity and stellar mass estimates. These systematic uncertainties are summarized in Table [3] 
First, we can mention that the stellar masses derived from photometric DR7 data can be systematic biased low due t o the offset 
between photometric and spectroscopic magnitudes. This bias can account for ~ 10% of uncertainty dFritz et alj|201 ll) . Then, we 
can highlight the impact of assuming one Schechter function to describe the overall distribution when one have a steep population 
faintward. This can affect the determination of M* and a parameters. When not using two functions, the determination of these 
parameters can systematic bias up the total luminosity determination by 6 - 15%, what will be directly translated to the M* 
estimates. Also related to the Schechter fit, we can adress the systematics of the covariance between M* and a parameters. If not 
taken into account, the uncertainties related to these parameters can biased low the total luminosity estimates by about 5%. The 
major systematic effect may be introd uced by the ado pt ed IMF when a ssuming a mass-to-ligh ratio to obtaine the stellar mass. 
For instance, a change from a standardlSarpeterj (jl955l) to lKroupal ( l2001l) IMF increses the M/Lby a factor of two. This translates 
into an increase by the same amonut in the stellar mass in our systems. Finally we can adress the ICL contribution. In this work, 
we did n ot co nsider the total light contribution to the stellar mass estimates. Thus, our values can be 2% systematic biased low 
(see Sect. ED. 



TABLE 3 

Main systematic uncertainties in stellar mass determination 



uncertainty 


% (approx.) 


photometric magnitude 


-10 


Scheschter fit 


+ 6-15 


covariance between M* and a 


-5 


IMF 


+100 


not considering ICL 


-2 




Fig. B.l. — Color-magnitude diagrams and luminosity function fits for our sample. 



